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I.  INTRODUCTION 


The  responses  of  a  millimeter-wave  (MMW)  radiometric  (passive)  sensor  and  a  MMW 
radar  (active)  sensor  to  varying  grassy  terrains  and  targets  were  measured  in  July  1980  by 
BRL  at  Aberdeen  Proving  Ground  (APG), Maryland.  These  measurements  were  made  as 
part  of  continuing  studies  of  MMW  sensors  for  smart  munitions.  The  analog  voltage 
signals  of  the  sensors  were  digitized  at  a  rate  of  2,000  observations  per  second, which  was 
the  same  rate  as  the  data  recorded  in  August  1978  at  Rome  ,  New  York.  Test  ranges  were 
selected  on  Spesutie  Island  at  APG.  Then  a  helicopter  was  instrumented  with  an  active 
and  passive  sensor  on  a  spinning  mount,which  revolved  at  three  (3)  revolutions  per 
second  (rps)  and  was  flown  over  the  ranges  on  runs  of  several  seconds  duration.  The 
passive  or  radiometric  sensor  was  similar  to  the  sensor  used  in  the  August  1978  teat  and 
thus  the  results  were  expected  to  have  similar  characteristics  (see  Figure  1 ). 

The  analysis  of  passive  and  active  sensors  is  an  on-going  investigation  at  BRL.  By  using 
the  Box-Jenkins1  time  series  modeling  approach  a  successful  simulation2  of  the  response 
of  the  passive  sensor  to  varying  grassy  terrains  in  the  1978  tests  at  Rome  has  provided  a 
realistic  way  to  introduce  false  alarms  and  the  distortion  of  the  target  signatures  due  to 
clutter  into  the  study  of  the  performance  of  various  weapon  systems.  The  response  of  the 
passive  sensor  was  modeled  by  an  Autoregressive-Intergrated-Moving-Average  (ARIMA) 
process,  which  was  first  order  autoregressive  and  first  order  moving  average.  The  purpose 
of  this  research  was  to  provide  a  simulation  of  the  response  of  the  passive  and  active 
sensors  used  in  the  1980  APG  tests  by  time  series  analysis.  The  model  for  the  passive 
sensor  was  expected  to  be  the  same  as  the  previous  study. 

On  the  other  hand,  the  signal  response  for  the  radar  (active)  sensor  was  different  in 
character,  in  that  the  signal  patterns  were  not  symmetric  about  the  meaa  Thai  is,  the 
signals  in  the  positive  directions  were  larger  than  in  the  negative  direction  (see  Figure  2). 
In  order  to  best  utilize  the  Box-Jenkins  time  series  modeling  approach,  a  data 
transformation  was  necessary.  The  purpose  of  the  data  transformation  was  to  employ  the 
data  in  a  form  that  will  improve  the  analysis,  without  altering  the  data  structure  (in  terms 
of  the  autocorrelation  and  partial  autocorrelation  function). 

One  of  the  Biomedical  computer  packages  (BMDP)  from  UCLA,viz.,  Statistical 
Software  for  Time  Series  Analysis,  was  used  for  this  analysis3.  There  are  two  estimation 
methods  of  the  model  parameters,  the  first  being  the  conditional  least  square  method 
(LSM)  and  the  second  being  the  backcasting  method  (BM). 


1  Box,G.  E  P.  and  Jenkins,  G.  M.,  "Time  Series  Analysis: 

Forecasting  and  Contror,Holden-Day,Inc.,San  Francisco,Ca.,1970. 

2  Maruyama,R.  T.,"A  Time  Series  and  Intervention  Analysis  of  Millimeter-Wave  (MMW) 
Radiometric  Data", ARBRL-TR-02338,  July  1981. 

3  Lon-Mu  Liu,"User’s  Manual  for  BMDQ2T(TSPACK)  Time  Series  Analysis  (Box- 
Jenkins)",Technical  Report  No.  57,  Department  of  Biomathematics ,  UCLA,  1979. 
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Plotted  Against  the 


.5  ms) 


IL  TIME  SERIES  ANALYSIS 


In  this  report  time  series  analysis  is  used  to  simulate  the  signals  from  a  passive  and  an 
active  (MMW)  sensor.  It  seems  fitting  to  explain  the  basic  ideas  of  time  series,  and  to 
introduce  some  terminology  so  that  the  reader  can  better  follow  our  report.  The  digitized 
values  of  the  analog  voltage  signals  represent  a  time  series  ZlvZ2,...,Z^  of  N  successive 
observations,  which  are  regarded  as  a  sample  realization  from  an  infinite  population  of 
such  time  series  that  could  have  been  generated  by  the  stochastic  process.  Both  the 
radiometer  and  radar  data  have  been  modeled  as  though  they  were  derived  from  a 
stationary  process,  which  implies  that  the  mean  and  variance  are  constant,  and  that  the 
error  terms  are  normally  distributed. 

An  autoregressive  model  is  one  in  which  the  current  value  of  the  process  is  expressed 
as  a  finite,  linear  aggregrate  of  the  previous  values  of  the  process  and  an  error  term  or  so- 
called  shock  a,  =  N(0,o-a2),  where  crfl  is  the  standard  deviation  of  the  white  noise.  If  the 
data  is  centered  about  the  mean  /i,  then  Z,  «•  Z,  —  /x  can  be  represented  as  an 
autoregressive  (AR)  process  of  order  p,  viz., 

2,  -  *,Z,_,  + </> 2^,-2  +  •  •  +  4>p2t-p  +  at  • 

The  model  contains  p+2  unknown  parameters  n,4> |,02*— »d>p.  <ra2  *  which  must  be 
estimated  from  the  data. 

An  autoregressive  operator  of  order  p  can  also  be  defined  in  terms  of  the  backward 
shift  operator  B  (  where  BZt  —  Zf_i  )  by 

4 '>(B)am  I— 4>2B2—  •  •  —$pBp  ,  (II.2) 

and  then  Equation  (II.  1 )  becomes 


4>{B)2t  -  a,  . 


(II.3) 


A  moving  average  (MA)  process  of  order  q  is  one  in  which  2,  is  made  to  be  linearly 
dependent  on  a  finite  number  q  of  previous  a ’s,  viz.. 


2,  -  a,  &iat— t  2  ‘  “  ^qat~v  *  (11-4) 

It  contains  q+2  unknown  parameters  /x,0lt...,09,  <r02,  which  must  be  estimated  from 
the  data. 

A  moving  average  operator  of  order  q  can  be  defined  by 

9{B)  -  1  -9yB  -  92B2  ->•••-  99B<  ,  (II.5) 

so  that  Equation  (II.4)  can  be  written  as 
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2,  -0(5)fl(  . 


(II.6) 


Frequently,  the  two  processes  are  combined  to  include  both  autoregressive  and  moving 
average  terms  in  the  model.  Thus,  the  mixed  autoregressive  -moving  average  model  is 

2,  —  0, 2t_x  +  •  •  •  +  0p2,_p  +  a,  —  0,a,_i  —  •  •  •  —  .  (H.7) 

or 


0(5)2,  -0(5)a,  ,  (II.8) 

where  there  are  p+q+2  unknown  parameters  #t;0 ),...,0p;  0i,...,0y;  <ra2  that  must  be 
estimated  from  the  data. 

Many  time  series  exhibit  nonstationary  behavior.  Certain  types  of  nonstationarity  can 
be  removed  by  differencing  the  data.  In  terms  of  the  backward  shift  operator  B,  the  time 
series  observations  (Z,)  can  be  differenced  d  times.  Thus,  the  new  operator  4>(5)  can  be 
defined  as 


4>(B)-0(B)(l-5)rf  ,  (II.9) 

where  0(5)  is  the  stationary  operator  in  Equation  (II.2).  The  general  model,  which 
represents  the  homogeneous  nonstationary  behavior  can  now  be  written  as 

<b(5)Z,  -  0(5)(l-5)rfZ,  -  0(5  )af  ,  (11.10) 

that  is 


0(5)^,  -0(5)a,  ,  (11.11) 

where 

W,  -(ISYZ,  .  (11.12) 

It  can  be  seen  that  Equation  (11.11)  means  that  the  d’th  difference  of  the  process  is 
stationary.  In  our  report  d  was  always  zero  since  our  original  data  was  stationary.  The 
process  defined  by  Equation  (11.11)  and  (11.12)  is  called  an  autoregressive  integrated 
moving  average  (ARIMA)  process  of  order  (p,d,q). 

Use  was  also  made  of  the  BMDQ2T(TSP ACK)  TIME  SERIES  PACKAGE  computer 
program  mentioned  in  reference  3.  This  program  evaluates  the  mean,  standard  deviation, 
autocorrelation  function  (ACF),  and  the  partial  autocorrelation  function  (PACT)  of  the 
N  successive  observations  of  the  Z,'s.  A  histogram  of  the  variable  Z  and  a  cumulative 
histogram  of  the  variable  Z  is  also  provided.  When  the  radar  data  was  analyzed  the  original 
data  (Z, )  was  transformed  into  (  Y,  )  by  a  natural  log  transformation  (Box-Cox)  and  then 
all  the  above  quantities  were  provided  for  the  (T, ). 
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The  ACF  is  a  diagnostic  tool,  which  tells  how  the  values  of  the  time  series,  (Z,  ZlHk), 
separated  by  k  time  steps  or  k  lags,  are  correlated  with  one  another.  The  PACF  is  a 
diagnostic  tool,  which  shows  the  autocorrelation  between  the  values  of  the  time  series, 
(Z,  Z,_*),  when  the  effects  of  the  variables  Z,_,  Z,_^  ....  ZlHk+1  are  fixed.  If  the  ACF 
and  PACF  of  the  original  data  are  plotted,  a  clue  is  given  for  a  tentative  starting  ARIMA 
model.  An  AR  process  of  order  p  has  an  ACF,  which  is  infinite  in  extent  with  damped 
exponentials  and/or  damped  sine  waves,  while  its  PACF  is  finite  in  extent  and  cuts  off 
after  p  nonzero  values.  An  MA  process  of  order  q  has  an  ACF,  which  is  finite  in  extent 
and  cuts  off  after  q  nonzero  values,  while  its  PACF  is  infinite  in  extent  and  dominated  by 
damped  exponentials  and/or  sine  waves.  A  mixed  process  of  order  (p,0,q)  has  an  ACF, 
which  is  infinite  in  extent  with  damped  exponentials  and/or  damped  sine  waves  after  the 
first  q-p  lags,  while  its  PACF  is  infinite  in  extent  and  dominated  by  damped  exponentials 
and/or  sine  waves  after  the  first  p-q  lags. 

After  a  tentative  model  has  been  identified,  the  ACF  and  PACF  of  the  residuals  (a, ), 
which  are  the  estimates  of  the  error  terms  for  a  particular  sampling  of  data,  are  plotted  to 
see  if  they  indicate  the  normality  or  white  noise  assumed  by  the  time  series  approach.  The 
residuals  of  a  white  noise  process  have  an  ACF  and  PACF,  which  have  value  1  at  zero  lag 
and  value  zero  elsewhere.  If  the  residuals  do  not  appear  to  be  normal,  a  new  ARIMA 
(p,d,q)  model  is  tried  until  the  residuals  indicate  normality.  The  program  also  estimates 
the  constants  in  the  (p,d,q)  model  by  the  conditional  least  squares  method  and  the 
backcasting  method.  For  our  problem  the  parameters  found  from  the  backcasting  method 
were  selected.  There  is  only  a  slight  difference  in  the  values  of  the  parameters  evaluated 
by  the  different  methods.  The  interested  reader  can  find  much  more  detail  on  time  series 
analysis  in  reference  1. 


Ill  MILLIMETER-WAVE  (MMW)  RADIOMETRIC  ANALYSIS 

A  previous  set  of  data,  which  was  taken  over  a  grassy  field  under  similar  conditions  at 
Rome,  N.Y.,  in  August  1978,  was  analyzed  in  1981  .  The  spin  rate  of  the  radiometer  at 
the  Rome  tests  was  4  rps.  The  result  of  that  investigation  was  that  an  ARIMA(  1,0,1) 
process  best  described  the  passive  signals.  Because  of  the  similarities  in  the  sensors  it  was 
expected  that  a  time  series  of  the  APG  data  would  again  result  in  a  first  order 
autoregressive-moving  average  (ARIMA(l,0,t)>  process. 

This  analysis  was  done  to  characterize  clutter  data  (background)  free  of  any 
intervention  such  as,  bodies  of  water, metallic  vehicles,metal  towers, etc.  Both  passive  and 
active  millimeter-wave  signals  were  collected  simultaneously.  The  sensors  were  mounted 
together  such  that  the  same  section  of  terrain  was  being  observed  at  the  same  instant  of 
time.  There  is  some  difference  in  the  area  of  the  ground  being  observed  by  the  two 
sensors.  The  footprint  of  the  active  sensor  is  smaller  than  that  of  the  passive  sensor. 

The  radiometric  data  shown  in  Figure  1  has  been  calibrated  in  Kelvins.  The  data  that 
the  following  analysis  concerns  are  uncalibrated.  To  obtain  calibrated  data  the  mean  (6.43) 
of  the  uncalibrated  data  must  be  substracted  from  each  uncalibrated  data  point  and  the 
result  multiplied  by  a  calbration  constant  of  81.4  K  per  unit. 
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The  Autocorrelation  Function  (ACF)  and  Partial  Autocorrelation  Function  (PACF)  of 
all  the  data  in  this  report  were  carried  out  to  36  lags  on  the  BMDP  package,  but  usually 
fewer  tags  were  needed  to  display  the  important  features.  For  the  passive  sensor  the  ACF 
and  PACF  are  shown  in  Figure  3,  and  the  values  plotted  are  given  in  Table  1.  The  decaying 
characteristics  of  the  ACF  are  theoretically  associated  with  an  AR(1)  model.  Thus,  an 
ARIMA(  1,0,0)  model  was  first  used  to  characterize  this  data. 

The  ARIMA(1,0,0)  fit  resulted  in  the  residuals  (a,)  of  this  fit  having  some  remaining 
structure  as  shown  in  Figure  4.  The  values  plotted  are  listed  in  Table  2.  The  ACF  at  lags 
one  (1)  and  two  (2)  as  well  as  the  PACF  at  lags  one  and  two  have  values  that  do  not 
indicate  normality.  The  lack  of  fit,determined  from  the  analysis  on  the  residuals,  indicates 
that  either  a  moving-average  term  or  an  autoregressive  term  should  be  added.  The  cutoff 
after  the  second  PACF  lag  indicates  that  an  autoregressive  term  should  be  added;  hence 
an  ARIMA(2,0,0)  was  attempted. 

The  ARIMA(2,0,0)  model  was  fitted  to  the  data  set  as  shown  in  Figure  5.  The  values 
plotted  are  given  in  Table  3.  The  residual  analysis  of  the  ARIMA(2,0,0)  model  did 
remove  the  significant  value  of  the  ACF  at  lag  one,  but  both  values  of  the  ACF  and 
PACF  at  lag  two  (2)  remain  significant,  indicating  a  need  for  further  analysis.  A  histogram 
of  the  residuals  in  this  case  is  shown  in  Figure  6,  and  gives  questionable  evidence  of  being 
normally  distributed.  The  only  unresolved  problem  area  in  the  analysis  of  the 
ARIMA(2,0,0)  was  with  the  remaining  residual  autocorrelation  at  the  second  lag. 

To  resolve  the  problem  of  the  significant  value  of  the  ACF  at  lag  2,  a  moving  average 
term  of  order  two  was  added  to  the  ARIMA(2,0,0)  structure.  The  AR1MA(2,0,(2))  fit, 
shown  in  Figure  7,  did  remove  both  the  large  values  at  lag  2  of  the  ACF  and  PACF  to  a 
point  where  it  was  adequate  for  the  purpose  of  characterizing  the  data.  The  values  are 
given  in  Table  4.  Therefore,  this  time  series  analysis  seemed  satisfactory. 

On  the  other  hand,  if  a  moving  average  term  were  added  to  the  ARIMA(  1,0,0)  model 
instead  of  another  autoregressive  term,  would  the  ARIMA(1,0,1)  structure  remove  the 
remaining  ACF  lags?  To  answer  this  question  the  the  ACF  and  PACF  of  the  residuals 
were  plotted  in  Figure  8,  and  their  values  are  given  in  Table  5.  A  histogram  of  the 
residuals  indicated  normality  as  can  be  seen  in  Figure  9.  The  only  negative  indication  is 
that  the  ACF  at  lag  two  (2)  and  the  PACF  past  lag  two  indicated  a  decaying  effect. 

The  analysis  of  these  residuals  indicated  that  an  additional  moving  average  term  might 
be  included  to  remove  the  autocorrelation  (ACF)  at  lag  two;hence,  the  ARIMA(  1,0,2) 
model  was  the  next  logical  step  in  the  sequence  of  analysis  for  the  millimeter-wave 
radiometric  data.  The  plot  of  the  residuals  shown  in  Figure  10  again  indicated  some  minor 
problems  with  the  analysis.  The  values  plotted  are  given  in  Table  6.  A  histogram  of  the 
residuals  seemed  normal  as  can  be  seen  in  Figure  11.  The  problem  is  that  the  value  of 
both  the  ACF  and  PACF  at  lag  three  is  now  believed  to  be  significant  Some  minor 
improvements,  such  as  a  decrease  in  the  standard  deviation,  have  taken  place  by  adding 
the  second  moving  average  term  to  the  ARIMA(1,0,1)  model,  but  these  improvements 
are  not  significant. 
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TABLE  1 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCOR RELATION 
FUNCTIONS  OF  RADIOMETRIC  DATA 


MEAN  -  6.43 

STANDARD  DEVIATION  -  .0595 


ACF 


670  Observations 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

r*  Lags  1-10 

.909 

.797 

.739 

.716 

.707 

.703 

.697 

.694 

.698 

.691 

11-20 

.680 

.673 

.668 

.656 

.643 

.628 

.606 

.590 

.583 

.577 

21-30 

.558 

.544 

.545 

.551 

.558 

.563 

.557 

.544 

.542 

.537 

31-33 

.519 

.499 

.494 

PACF 

670  Observations 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0*  Lags  1-10 

.909 

-.177 

.281 

.086 

.127 

.100 

.066 

.097 

.103 

-.010 

11-20 

.070 

.040 

.032 

-.010 

.017 

-.027 

-.046 

.012 

.007 

-.017 

21-30 

-.074 

.038 

.044 

.019 

.063 

.038 

-.007 

.000 

.080 

-.024 

31-33 

-.049 

-.019 

.055 
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TABLE 2 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RESIDUALS 
(ARIMA(1,0,0)) 


MEAN  -  -.0000301 
STANDARD  DEVIATION  -  .0245 


ACF 

670  Observations 


1 

2 

3 

4 

5 

6 

7 

8 

9  1 

!  n  Lags  1-9 

! 

i 

.169 

-.262 

-.159 

-.044 

.006 

.043 

.019 

-.019 

.095  \ 

i 

i 

i 

i 

i 

i 

i 

PACF 

'  670  Observations 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

<f»k  Lags  1-9 

.169 

-.299 

-.057 

-.087 

-.034 

.006 

-.015 

-.016 

.121 

j  i 


i  t 

w 
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TABLE  3 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RESIDUALS 
(ARIMA(2,0,0)) 


MEAN  -  -.0000174 
STANDARD  DEVIATION  -  .0241 


ACF 

670  Observations 

12  3456789 


rk  Lags  1-9  .053  -.242  -.070  .016  .035  .068  .042  -.019  .116 


PACF 

670  Observations 

12  3456789 


<f>k  Lags  1-9  .053  -.246  -.044  -.040  .009  .061  .047  .011  .157 
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TABLE 4 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RESIDUALS 
(ARIMA(2,0,(2))) 


MEAN  -  -.0000403 
STANDARD  DEVIATION  -  .0229 


670  Observations 


ACF 

234S6789 


rk  Lags  1-9  -.005  .020  -.172  -.041  -.068  -.008  .004  -.052  .086 


PA  CF 


670  Observations 

123456789 


<f>k  Lags  1-9  -.005  .020  -.172  -.043  -.064  -.038  -.009  -.079  .074 
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TABLE  5 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RESIDUALS 
(ARIMA(  1,0,0) 


MEAN  -  -.0000150 
STANDARD  DEVIATION  -  .0237 


ACF 


670  Observations 

_ 1  2  3  4  5  6  7  8  9 


rk  Lags  1-9  -.047  -.141  -.045  .021  .036  .066  .059  -.023  .128 


PACF 


670  Observations 

_ 1  23456789 

<t>k  Lags  1-9  -.047  -.144  -.061  -.005  .022  .071  .079  .009  .158 
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Figure  9.  Histogram  of  Residuals  of  the  ARIMA(1,0,1 ) . 


TABLE 6 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RESIDUALS 
(ARIMA(  1,0,2)) 


MEAN  -  -.0000429 
STANDARD  DEVIATION  -  .0231 


ACF 

670  Observations 

_ 123456789 

rk  Lags  1-9  .087  .025  -.194  -.068  -.092  -.022  -.013  -.047  .075 


PACF 

670  Observations 

123456789 


<f>k  Lags  1-9  .087  .018  -.199  -.036  -.075  -.046  -.024  -.082  .068 


29 


Nevertheless,  the  ARIMA(1,0,3)  model  was  tentatively  used  to  remove  the  significant 
values  of  both  the  ACF  and  PACF  at  the  third  lag.  A  plot  of  the  ACF  and  PACF  of  the 
residuals  is  shown  in  Figure  12,  with  the  values  given  in  Table  7.  After  all  these  modeling 
attempts,  the  ARIMA(  1,0,3)  model  seems  to  characterize  the  radiometric  data  best.  Both 
the  ACF  and  PACF  of  the  residuals  seem  small  and  the  standard  deviation  of  0.0224.  is 
smaller  than  any  other  model  that  was  studied  as  can  be  seen  in  Table  8. 

At  this  point  in  the  analysis  a  reflection  on  the  ability  of  each  tentative  model  to 
characterize  this  data  parsimoniously  should  be  addressed.  Although  the  ARIMa(  1,0,3) 
may  be  slightly  better  than  the  other  models,  a  summary  of  the  ARIMA  models  in  terms 
of  parameters  and  residuals  indicates  that  the  ARIMA(1,0,1)  structure  is  almost  as  good 
as  any  other  satisfactory  model,  and  it  only  requires  two  parameters.  Therefore,based  on 
the  principle  of  parsimony  the  ARIMA(1,Q,1)  model  was  chosen  to  characterize  the 
radiometric  data. 

If  the  following  definitions  are  made 

a,  ~N(0,(ru2),  (111.1) 

and 

2, -Z, -M,  (III.2) 

where  the  mean  of  Z,  is  n  ,  then  the  ARIMA(  1,0,1 )  becomes 

(  1  -  4»XB  )2,  -  (  1  -  9}B  )at  ,  (111.3) 

where  the  backward  shift  operator  B  is  defined  as  BZ,  —  Z,_,  ,  and  Ba,  —  a,_,  .  If 
Equation  (III.3)  is  expanded  and  rearranged,  then  in  terms  of  the  original  data  Z,  becomes 

Z,  -  (  1  -  0,  )m  +  *,Z,_,  +  a,  —  0, a,_,  .  (II  1.4) 


IV  MILLIMETER-WAVE  (MMW)  RADAR  ANALYSIS 

The  radar  data  was  collected  simultaneously  with  the  radiometer  data  at  APG,  Md.  The 
radar  is  an  active  sensor,  which  means  that  it  sends  out  a  signal  that  is  scattered  otY  objects 
on  the  ground  and  back  into  its  receiving  antenna.  On  the  other  hand  the  radiometer  does 
not  send  signals  ,  but  only  receives  them  from  the  ground;  hence  it  is  referred  to  as  a 
passive  sensor.  Uncalibrated  signals  from  the  radar  sensor  are  presented  in  Figure  1 3.  To 
obtain  the  calibrated  radar  data  of  Figure  2  the  zero  cross  section  reference  level  (2.95)  of 
the  uncalibrated  data  must  be  substracted  from  each  uncalibrated  data  point  and  the  result 
multiplied  by  a  calibration  constant  of  8.0  square  meters  per  unit. 

Figure  14  is  a  histogram  of  the  uncalibrated  radar  data,  which  shows  the  skewness  of 
these  signals.  It  is  obvious  that  the  deviations  are  not  symmetric  about  the  mean  as  well  as 
having  positive  deviations  larger  in  magnitude  than  those  in  the  negative  direction. 
Because  of  this  skewness  in  the  radar  data,  a  data  transformation  was  necessary  to  best 
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TABLE 7 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RESIDUALS 
(ARIMA(1,0,3)) 


MEAN  -  -.0000469 
STANDARD  DEVIATION  -  .0224 


670  Observations 


ACF 

234  5  6789 


rk  Lags  1-9  .012  .025  -.001  -.067  -.018  -.017  .010  -.063  .079 


PACF 


670  Observations 

1234  5  6789 


4>k  Lags  1-9  .012  .025  -.001  -.068  -.017  -.013  .011  -.068  .078 
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TABLE  8 


SUMMARY  OF  ARIMA(*)  MODELS  ENTERTAINED  FOR  PASSIVE  SENSORS 

(BACKCASTTNG  METHOD) 


ARIMA(p,d,q) 

(p.d,q) 

ESTIMATED 

PARAMETERS 

WHITE  NOISE 

NOTE 

(1) 

(1.0.0 

0,  -0.9127 

a - .00003 

Some  remaining  structure  in  the 

centered 

Sa+  —  0.0245 

residuals  [AR(2)  or  MA(1)]. 

(2) 

(2,0,0) 

0,  - 1.086 

a  --.000017 

ACF  at  Lag  2. 

centered 

02 - -1902 

S„  -  0.0241 

(3) 

(2.0.C2)) 

0,  -  1.0913 

a  ——.00004 

ACFat  Lag  3. 

centered 

£2-  -1274 

Sa  -  0.0229 

02  -  0.3866 

(4)* 

(1.0.1) 

0,  -  0.8537 

a- -.00001 5 

ACF  at  Lag  2  with  PACF 

centered 

fl,  -  -.3569 

Sa  -  0.0237 

at  Lag  2  with  decay. 

(5) 

(1.0,2) 

0,-0.9676 

S  -  -.000043 

ACFat  Lag  3. 

centered 

0, - .028 

Sa  -  0.0231 

02  -  0.4067 

(6) 

(1.0.3) 

0,  -  0.9874 

a  -  -.000947 

centered 

0,  -  -.0654 

Sa  -  0.0224 

?2- 0.4342 

0.2492 

+ The  symbol  Sa  is  the  estimated  standard  deviation. 

*The  ARTMA(  1,0,1)  results  are  as  good  as  any  above,  and  based  on  the  principle  of  parsi¬ 
mony  one  would  choose  the  p— 1,  d— 0,  q— 1  model. 
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Figure  13.  The  Actual 
is  Plotted  A 


Figure  14.  Histogram  of  the  Actual  Uncalibrated  Radar  Data 


utilize  the  Box  and  Jenkins  method  on  this  data  set. 

A  natural  logarithm  (In)  transformation  was  selected  to  transform  the  radar  data  (Zt) 
into  a  new  variable  O',)  in  such  a  way  that  the  new  deviations  are  symmetric  abour  their 
mean.  This  logarithm  transformation  maintains  the  order  of  the  original  data  and 
symmetrically  transforms  the  data  to  enhance  the  time  series  analysis. 

The  data  transformation  used  is  based  upon  a  simple  idea  of  equating  the  maximum 
and  minimum  distance  from  the  mean  value.  Hence,  take  the  maximum  logarithm  value 
of  the  observations  minus  the  mean  value  and  equate  it  to  the  mean  logarithm  value 
minus  the  minimum  value.  This  procedure  equates  the  magnitude  of  the  maximum  and 
minimum  deviations  from  the  mean,  thereby  improving  the  symmetric  properties. 


max  //»(•)  —  mean  Ini •)  —  mean  /«(•)  —  min  //»(•)  .  (IV.l) 

or 

Ini  3.40  +  A  )  -  Ini  3.08  +  A  )  -  Ini.  3.08  +  A  )  -  Ini  2.96  +  A  )  .  (IV.2) 

The  solution  for  A  is  -2.89. 

The  transformation,  Y,  —  Ini  Z,  —  2.89  ),  which  is  based  on  improving  the  symmetric 
properties,  needs  a  second  minor  adjustment.  This  minor  adjustment  was  a  correction  in 
the  value  of  the  parameter  A  to  -2.91. 

Yt-lniZt-  2.91),  (IV.3) 

where 


mean  Y,  —  —1.85  , 


and 


maximum  Yt  —  —.71 34  , 
minimum  Y,  —  —2.995  , 
standard  deviation  SY  —  0.370  . 


More  details  of  the  evaluation  of  the  parameter  A  are  given  in  Appendix  A. 

A  nice  property  of  this  ln(  •  )  transformation  is  that  the  order  of  the  observations 
remain  the  same  as  the  original  data.  Hence,  the  ACF  and  PACF  of  the  original 
observation  as  well  as  the  transformed  data  are  similar.  This  similarity  can  be  seen  by 
comparing  Figure  15  with  Figure  16.  The  values  plotted  are  listed  in  Table  9  and  Table 
10,  respectively. 
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TABLE  9 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RADAR  DATA 


MEAN  -  3.08 

STANDARD  DEVIATION  -  .0697 


ACF 


670  Observations 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

rk  Lags  1-10 

.898 

.682 

.457 

.289 

.189 

.135 

.113 

.111 

.121 

.140 

11-20 

.162 

.187 

.213 

.239 

.256 

.252 

.220 

.178 

.141 

.127 

21-30 

.134 

.150 

.159 

.152 

.138 

.130 

.135 

.143 

.153 

.162 

31-33 

.164 

.147 

.112 

PACE 

670  Observations 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

4>k  Lags  1-10 

.898 

-.647 

.253 

.081 

.062 

.009 

.082 

016 

.028 

.054 

11-20 

.045 

.057 

.054 

.048 

.001 

-.046 

-.019 

.081 

-.019 

.097 

21-30 

-.011 

.002 

-.003 

-.055 

.101 

-.004 

.034 

-.046 

.096 

-.019 

31-33 

-.040 

-.044 

-.013 
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TABLE  10 


ESTIMATED  AUTOCORRELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  TRANSFORMED  RADAR  DATA 


MEAN  ”  -1.85 

STANDARD  DEVIATION  -  .370 


ACF 


670  Observations 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

rk  Lags  1-10 

.895 

.689 

.481 

.330 

.242 

.196 

.178 

.171 

.168 

.171 

11-20 

.180 

.194 

.215 

.237 

.250 

237 

.198 

.157 

.125 

.117 

21-30 

.132 

.158 

.179 

.179 

.171 

.162 

.158 

.154 

.153 

.154 

31-33 

.151 

.137 

.111 

PACF 

670  Observations 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

<f>k  Lags  1-10 

.895 

-.567 

.164 

.103 

-.018 

.022 

.054 

-.002 

.029 

.079 

11-20 

.002 

.065 

.071 

.022 

-.006 

-.070 

.006 

.060 

-.030 

.096 

21-30 

.025 

.014 

-.007 

-.041 

094 

-.026 

.025 

-.024 

.054 

-.002 

31-33 

-.009 

-.028 

-.035 
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The  sinusoidal  pattern  in  the  ACF  and  the  cutoff  pattern  in  the  PACF  indicate  some 
type  of  autoregressive  (AR)  process.  A  second  order  autoregressive  model  was  therefore 
entertained.  The  values  of  the  ACF  and  PACF  for  the  residuals  are  plotted  in  Figure  17, 
and  displayed  in  Table  11. 

Although  Table  1 1  indicates  a  reasonable  degree  of  success,  the  residuals  plotted  as  a 
histogram  in  Figure  18  indicated  that  the  residuals  were  not  gaussian  (normal).  Hence,  the 
residuals  did  not  satisfy  the  basic  assumption  in  the  Box  and  Jenkins  method  that  is 
explained  in  more  detail  in  their  book1. 

To  improve  upon  the  time  series  model  a  number  of  moving  average  terms  were 
introduced  into  the  analysis  to  remove  any  significant  value  of  the  ACF  and  PaCF  that 
may  exist  at  the  different  lags,  as  well  as  to  satisfy  the  underlying  white  noise  assumptions 
of  normality  (N(0,<ra2)).  See  Table  12  for  a  breakdown  of  the  ARIMA  time  series 
analysis  attempted  during  this  investigation. 

The  ARIMA(2,0,(2,3,4))  process  was  able  to  characterize  the  millimeter  wave  radar 
data  well,  where  the  random  shock  was  truly  random  and  symmetric  as  Figure  19 
indicates.  The  question  that  always  arises  in  this  type  of  modeling  effort  is  how  general  is 
this  particular  ARIMA  model.  In  order  to  answer  such  a  question  ,  some  sets  of 
observations  other  than  those  already  reported  on  were  analyzed  with  the  same  ARIMA 
model.  The  results  are  reported  in  Table  1 3. 

This  particular  time  series  structure  was  able  to  characterize  the  other  data  sets,  where 
the  estimated  parameters  remained  relatively  the  same,  indicating  some  degree  of 
generality  of  this  time  series  process. 

If  the  following  definitions  are  made, 

a,  =Af(0,<7a2),  UV.4) 

and 

Wt  -  Y,  -  n  ,  (IV.5) 

where  the  mean  of  Y,  is  n ,  then  the  ARIMA(2,0,(2,3,4))  becomes 

(  1  -  -  <f>2B2  )W,  -  (  1  -  82B2  -  93B 3  -  9aBa  )a,  ,  (IV.6) 

where  the  backward  shift  operator  B  is  defined  as  BVY,  “  W,_x  ,and  Bo,  —  of_,  .  If 
Equation  (IV.6)  is  expanded  and  rearranged, 

Wt  -  <f>xW,_x  +  2  +  a,-  820,-2  ~  830,-3  -  840,-4 ,  (IV.7) 

thus,  Y,  becomes 

Y,  -  n  +  *,(  Y,-X -n)  +  4>2{Y,-2-n)  +  a,- 9^,-2  ~  830,-3  -  8 40,^  .  (IV.8) 

Finally,  in  terms  of  the  original  data  Z,  becomes 

Z,  -eV'+2.91.  (IV .9) 
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ACF 


3 


TABLE  11 


ESTIMATED  AUTOCOR RELATION  AND  PARTIAL  AUTOCORRELATION 
FUNCTIONS  OF  RESIDUALS  OF  TRANSFORMED  RADAR  DATA 

(ARIMA(  2,0,0)) 


MEAN  =  .0000668 
STANDARD  DEVIATION  -  0.133 


ACF 


670  Observations 


1 

2 

3 

4 

5 

6 

_ 7 _ 

8 

9 

10 

rk  Lags  1-10 

.110 

-.030 

-.101 

.012 

.101 

.060 

.087 

.084 

.008 

.057 

11-20 

.044 

.015 

.022 

.016 

.119 

.112 

.008 

.048 

-.031 

.001 

21-24 

.007 

.031 

.161 

-.008 

PACF 

670  Observations 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

<J>k  Lags  1-10 

.110 

-.073 

-.088 

.030 

.087 

.033 

.094 

.091 

.006 

.075 

11-20 

.040 

-.006 

.016 

.003 

.096 

.081 

-.011 

.065 

-.038 

-.018 

21-24 

-.019 

-.006 

.121 

-.054 
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TABLE  12 


SUMMARY  OP  ARIMA  MODELS  FOR  TRANSFORMED  RADAR  DATA  BY  BACKCASTING  METHOD 

(Y,  - ln(Z ,  —  2.91);K  ~-1.85,Sy  -0.37) 


ARlMA(p,d,q) 

(p,d,q) 

ESTIMATED 

PARAMETERS 

WHITE  NOISE 

NOTE 

(1) 

(2.0.0) 

0,  -  1.4309 

a  -  .000067 

Significant  value  for  ACFand  PACF  at  Ljg(l). 

Centered 

02  -  -.5962 

Sa  -0.133 

(2) 

(2.0.1) 

0,  -  1.2816 

a - .000061 

Residuals*  not  normal,N(0,ora  2). 

Centered 

02  -  -.4638 

5^  —  0.131 

0,  “  -.2435 

(3) 

(2,0,2) 

0,-1.0331 

a  -  -.000061 

Residuals*  not  noimal,N(0,cro2) 

Centered 

02  -  -.2715 

Sa  -0.130 

0,  —  —.4995 

02  -  -.2381 

(4) 

(2,0,  (2,3)) 
Centered 

0,-1.5021 

02  --5787 

02 -.1790 

03  -  .2445 

a  -  -.000061 
Sa  —  .131 

Residuals*  not  norma],N(0,o-o2) 

Significant  value  of  ACF  and  PACF  at  Lag(4). 

(5) 

(2,0,  (2, 3,4)) 
Centered 

0,  - 1.5261 

02  --.5339 

02  -  .2920 

03  -  .3870 

04  -  .2007 

a  --.00112 

Sa  -0.129 

Residuals  are  normal,N(0,®-fl2). 

The  ACF  and  PACF  are  well  behaved. 

+BMDQ2T  (TSPACK)  Time  Series  Package;  Program  Revised  July  1,  1979  .Copyright  (c) 
1979,Regents  of  University  of  California.  (Storage  16000.) 

*The  histogram  of  the  residuals  indicates  that  the  normality  assumption  for  the  Box  and 
Jenkins  ARIMA  process  does  not  hold. 
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Figure  19.  Histogram  of  Residuals  of  ARIMAC2 , 0 , C2 , 3 ,4) > 

for  the  Transformed  Radar  Data. 


TABLE  13 


OTHER  DATA  SETS 
(BACKCASTI no  method) 


ARIMA(p,d,q) 

(p»d,q) 

ESTIMATED 

parameters 

WHITE  NOISE 

NOTE 

(1) 

(2,0,(2,3,4)) 

0,  ”  1.5261 

a-  -.00112 

Histogram  of  Residuals  “  N(Q,cra2)* 

Centered 

02  «  -.5339 

Sa  -0.129 

The  ACF  end  PACF  arc  -acI!  behaved  . 

Y  -  -1.85 

d2  -  0.2920 

Sy  -0.370 

flj  0.3870 

04  «  0.2007 

(2) 

(2,0,(2,3,4)) 

0i  —  1 .5273 

a  -  -.00121 

-N(0,cra2) 

Centered 

02  —  —.5360 

Sa  -0.126 

The  ACF  end  PACF  arc  well  behaved. 

Y  —  —1.88 

02  1:3  0.3807 

Sy  -0.338 

03  -  0.3505 

04  —  0.1469 

(3) 

(2,0,(2,3,4)) 

0,-1.5026 

a  -  -.00049 

-N(0,cra2) 

Centered 

02  -  -.5141 

Sa  -0.134 

The  ACF  and  PACF  are  well  behaved. 

Y  —  —1.80 

02  -  0.2972 

Sy  -0.404 

03  -  0.3682 

04  —  0.1162 

s - 

In  all  three  cases  the  histograms  of  the  residuals  of  the  ARIMA  fit  indicated 
normality. 
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V  SUMMARY 


The  purpose  of  this  research  on  the  1980  APG  test  data  was  to  provide  a  simulation  of 
the  response  of  the  passive  and  active  sensors  by  time  series  analysis.  Two  successful 
models  have  been  developed.  These  simulations  have  provided  a  realistic  method  to 
introduce  false  alarms  and  the  distortion  of  target  signatures  due  to  the  clutter  from  grassy 
fields  into  the  study  of  the  performance  of  various  weapon  systems.  The  important 
findings  of  this  research  are  presented  below. 

The  millimeter  wave  radiometric  data,  collected  in  July  1980  at  APG  and  analyzed  in 
this  report,  has  resulted  in  the  same  ARIMA(l,0,l)  structure  as  the  August  1978 
radiometric  data  collected  at  Rome,  New  York,  with  a  similar  sensor.  A  plot  of  (L ,Zf )  for 
L— 1,500  is  shown  in  Figure  20,  where  Z,  is  generated  using  Equation  (IIJ.4)  and  the 
parameter  values  in  Table  8,  and  then  calibrated  in  Kelvins.  If  this  plot  is  compared  with 
the  actual  radiometric  data  in  Figure  1,  the  agreement  is  good  except  for  the  seasonal 
variation  of  the  actual  data  with  a  period  of  approximately  300  lags. 

There  are  some  possible  physical  explanations  of  this  cyclic  behavior,  which  were  not 
present  in  the  Rome  tests,  viz., 

1.  In  the  Rome  tests  the  radiometer  was  the  only  sensor  being  tested, 
whereas  in  the  1980  APG  tests  the  radiometer  was  placed  on  a  spinning 
mount  along  with  an  active  sensor.  The  two  sensors  shared  the  same 
receiving  antenna,  and  there  may  have  been  some  radar  interference 
with  the  radiometer. 

2.  There  may  have  been  some  problem  with  the  slip  rings  through 
which  the  radiometric  signal  was  carried  to  the  recording 
devices. 

3.  There  may  have  been  some  radio  frequency  interference,  (RFI). 

To  address  the  cyclic  behavior  of  the  data  by  time  series  modeling  is  a  non  trivial 
problem  of  seasonality,  which  has  not  been  attempted  in  this  report  Basically,  our 
simulation  is  felt  to  be  an  adequate  representation  of  the  original  data. 

There  are  some  other  differences  between  the  Rome  tests  and  the  APG  tests.  Both 
data  sets  were  digitized  at  the  rate  of  2000  samples  per  second.  Since  the  spin  rate  of  the 
radiometer  was  4  rps  at  Rome  and  3  rps  at  APG,  there  were  500  observations  per 
revolution  of  the  sensor  at  Rome  compared  with  =  667  observations  per  revolution  at 
APG.  In  each  test  the  helicopter  flew  at  approximately  the  same  speed.  This  means  that 
the  consecutive  digitized  footprints  of  the  sensor  on  the  ground  were  closer  to  each  other 
at  APG  than  at  Rome.  Therefore,  the  ACF  of  the  APG  data  had  a  stronger  autoregressive 
structure;  i.e.,  the  values  of  the  ACF  at  the  various  lags  were  closer  to  unity  for  the  36 
tags  studied  than  the  values  of  the  ACF  for  the  Rome  data. 
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The  Simulated  Contrast  Temperature  Data  T(K)  for  a  Field  at  90m  Slant 
Range  is  Plotted  Against  the  Integer  Number  L  of  Time  Steps,  (A r  -  .Sms) 


Finally,  the  two  tests  were  not  carried  out  on  identical  grassy  fields,  nor  were  the 
weather  conditions  exactly  the  same.  Nevertheless,  it  is  worth  noting  that  the  same 
ARIMA(  1,0,1)  model  adequately  describes  both  the  Rome  and  APG  data. 

On  the  other  hand  the  millimeter  wave  radar  data,  collected  simultaneously  with  the 
July  1980  radiometer  data  at  APG  and  analyzed  in  this  report,  at  first  had  to  be 
transformed  by  a  Box-Cox  transformation  before  the  time  series  analysis  could  be  done. 
The  reason  for  the  transformation  was  that  the  the  original  radar  data  was  skewed  and  the 
deviations  were  not  symmetric  about  the  mean  as  well  as  having  positive  deviations  larger 
in  magnitude  than  those  in  the  negative  direction.  In  order  to  balance  the  positive  and 
negative  magnitudes  of  the  deviations,  the  properties  of  a  log  transformation  with  a  shift 
of  -2.91  was  used.  Then,  this  transformed  data  was  analyzed  using  the  Box  and  Jenkins 
process.  The  analysis  indicated  an  autoregressive-  moving  average  ARIMA  model  of  order 
p—2  and  q— 2,3,4  .  A  plot  of  (L £,)  for  L— 1,500  is  shown  in  Figure  21.  where  Z,  is 
generated  using  Equations  (IV.8)  and  (IV.9)  and  the  parameter  values  in  Table  12,  and 
then  calibrated  as  radar  cross  section  in  meters  squared.  If  this  plot  is  compared  with  the 
original  data  of  Figure  2,  the  agreement  is  quite  satisfactory. 
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Simulated  Radar  Cross  Section  trim1)  for  a  Field  at  90m  Slant  Range 
otted  Against  the  Integer  Number  L  of  Time  Steps,  (A t  —  .Sms)  . 
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APPENDIX  A 


A  time  series  data  transformation  that  is  required  for  the  radar  data  in  using  the  Box- 
Jenkins  modeling  approach  is  as  follows: 

Data  Transformation  (Box-Cox) 

•  Order/ranking  of  original  observations  is  maintained  after  the  transformation. 

*  The  deviations  of  the  transformed  data  are  symmetric  about  the  mean. 

Hence, a  natural  logarithmic  transformation  /«(•)  was  selected,  which  maintains  the 
order  of  the  original  observations.  The  minimum,  mean,  maximum,  and  standard 
deviation  of  the  random  variable  {Z, }  for  the  first  data  set  is: 

min  ln(Zt )  —  1.09  , 

mean  IniZ f )  -*  1.12  , 

max  In  (Z, )  —  1 .22  , 

standard  deviation  of  //»(•)  —  0.0223 

One  method  of  accomplishing  symmetry  about  the  mean  is  by  setting  the 
transformation  of  the  original  data  { Z,)  such  that  the  extreme  deviations  from  the  mean 
are  equal.  That  is, 

t  Ini  3.40  +  A  )  -  Ini  3.08  +  A  )  ]  -  [  Ini  3.08  +  A  )  -  Ini  2.96  +  A  )  ]  (A.1) 

where  A  is  a  shift.  Then 

Ini  3.40  +  A  )  +  /«(  2.96  +  A  )  —  2  //i(  3.08  4-  A  ) .  (A.2) 

If  the  //»(•)  is  removed, 

(  3.40  +  A  X  2.96  +  A  )-(  3.08  + A  )* ,  (A.3) 

where  A  “  —2.89  . 

Using  the  { Ini  Z,  —  2.89  )  }  transformation  results  in 

min  Yt  “  —2.66 , 
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I 


mean  Y,  ——1.73 , 
max  Y,  —  —0.673 , 
standard  deviation  —  0.329 


and 


A^  —  mean  Yt  —  min  Yt  —  0.93, 

A M  —  max  Yt  —  mean  Y,  —  1.057  . 


Then, 

~  “  -127* 

which  implies  that  the  deviation  in  the  maximum  (positive  direction)  is  larger  by  .127.  It  is 
desirable  to  back  shift  by  .5(.127)  —  .0635  .  So  an  adjustment  was  made  to  a  shift  of  -2.91 
with  the  result  that 


min  Y ,  —  2.9957  , 
mean  Y,  ——1.85 , 
max  Yt  ——0.71335 , 
standard  deviation  —  0.370 


and 


Al  -  1.145 , 
Aw  - 1.138 . 


Hence,  the  Y,  —ln(  Zt  —  2.91  )  data  transformation  accomplishes  the  objective  of  making 
the  deviation  symmetric  about  is  mean. 
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